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Abstract 

We explicitly take into account the effect of hydrodynamic expansion prohle on the 
gluonic breakup of J/'0’s produced in an equilibrating parton plasma. Attention is 
paid to the space-time inhomogeneities as well as Lorentz frames while deriving new 
expressions for the gluon number density n^, average dissociation rate (P), and "0 survival 
probability S. A novel type of partial wave interference mechanism is found to operate 
in the formula of (P). Nonrelativistic longitudinal expansion fro small length of the 
initial cylinder is found to push the S{pt) graph above the no flow case considered by us 
earlier [9]. However, relativistic flow corresponding to large length of the initial cylinder 
pushes the curve of S'(pt) downwards at LHC but upwards at RHIC. This mutually 
different effect on S'(pr) may be attributed to the different initial temperatures generated 
at LHC and RHIC. 


PACS numbers: 12.38M 
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I INTRODUCTION 


Extensive literature exists on the possible suppression [1]- [ 6 ] of the J/'ip mesons in a 
quark-gluon plasma and their proposed regeneration [7]. Among the well known mech¬ 
anisms of J/ip dissociation the one due to gluonic bombardment [ 8 ] deserves special 
attention here. Recently the present authors [9] considered the statistical mechanics 
of important physical observables viz. the gluon number density, thermally-averaged 
g — Ip break-up rate, and the ip meson survival probability appropriate to RHIC/LHC 
initial conditions. We found [9] that these observables are significantly affected if one 
employs improved expressions for the gluon distribution function, g — ip relative flux, 
and Ip meson formation time. 

Of course, it is a well-recognized fact that the longitudinal/transverse expansion 
of the medium controls the master rate equations [ 10 ] for the time-evolution of the 
plasma temperature and parton fugacities. But the literature does not tell how the fluid 
velocity prohle itself influences the Lorentz transformations connecting the rest frames 
of the hreball, plasma, and ip meson. In other words, since the flow velocity profile 
causes inhomogeneities in space-time, hence the scenario oi .7/ ip gluonic break-up may 
be affected in a quite nontrivial manner and the aim of the present paper is to address 
this hitherto unsolved problem. 

Sec.II below recalls a few aspects of relativistic hydrodynamics for the sake of ready 
reference. Next, a new expression for the gluon number density is derived in Sec.Ill 
showing how an extra dilation factor 7 associated with the flow appears. Next, in Sec.IV 
careful Lorentz transformations are used to calculate the flux-weighted cross section and 
explicit dependence is brought out on the hydrodynamic velocity w observed in 7/ip rest 
frame. Next, Sec.V develops the machinery for computing the 7/ip survival probability 
as a function of transverse momentum. Finally, Sec. VI summarizes our main conclusions 
applicable to nonrelativistic/relativistic flows. 


II ASPECTS OF HYDRODYNAMICS 

Preliminaries : Consider the rest frame of the hot, dense hreball produced in ultrarela- 
tivistic heavy ion collision. Within an initial time span = tq ~ 0.5 frm/c it is supposed 
to achieve local thermal equilibrium. The plasma now expands rapidly, gets cooled at 
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the expense of internal energy, and is driven towards chemical equilibration through par- 
tonic reactions. The plasma’s life ends, i.e., freeze-out occurs at the instant tufe when the 
temperature drops to 200 MeV, say. Such a picture of collective flow is known to have 
profound effect on the measured dilepton spectrum [ 11 ], dependence of colour screening 
mechanism on equation of state [ 12 ], anisotropy in the transverse-momentum distri¬ 
bution [13] of output hadrons, azimuthal asymmetry of J/ijj suppression in noncentral 
heavy-ion collisions [14], etc. 

Equation of motion; We employ the units h = c = 1 and follow closely the hydro- 
dynamic summary given by Pal et al [12] based upon cylindrical symmetry appropriate 
to central collisions. In the fireball frame a general time-space point x and the fluid 4 
velocity u have the form 

x = {t,x) ] u = (7,7^) 

7 = ( 1 + 5 ")‘'" = ( 1 -« 7 '‘’' ( 1 ) 


where v is the local 3 velocity and 7 is the corresponding Lorentz factor. Ignoring 
viscosity the conservation law for the energy- momentum tensor reads 


= 0 ; = (e + 


( 2 ) 


where the energy density e and pressure P are supposed to be measured in a frame 
comoving with the plasma. The relationship between the fireball usual time t and medium 
proper time r is, of course 


dr 

dt 


1 

) ti ^ t ^ hife 

7 


( 3 ) 


Longitudinal Expansion; In Bjorken’s boost-invariant, one-dimensional flow the pro- 
hle admits a simple analytical solution 

v = ; T = {t^ - > To (4) 

where is a unit vector along the collision axis. The corresponding temperature T oc 
is known to fall rather slowly, e.g. in the case of gold nuclei colliding at RHIC. 
Hence the J/fi suppression remains comparatively high viaihe colour screening/gluonic 
dissociation mechanisms. 

Transverse Expansion; As summarized in Refs. [12]-[11] the 4-velocity prohle be¬ 
comes quite intricate and a numerical integration of the dynamical equations (2,3) be¬ 
comes very hard in the case of 3 -|- 1 dimensional expansion. For our purpose it will 
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suffice to assume that the collective flow occurs only along the lateral directions without 
rotation as given by the empirical ansatz 

.2n1/2 


V = -Cr 

t 


; r={t^-ry 


( 5 ) 


where ir is the unit vector along r in a cylindrical coordinate system (r, z,(l)). It is known 
from numerical solution of ( 2 ) that transverse flow causes very rapid cooling in the case 
of lead nnclei colliding at LHC. Hence the J/'ip survival chance becomes relatively high 
niathe colonr screening/glnonic dissociation scenarios. We now proceed to formnlate the 
statistical mechanics of some physical observables following closely the logic of Ref. [9]. 


Ill GLUON NUMBER DENSITY 


Preliminaries : Working in the fireball rest frame and assnming local thermal eqnilib- 
rinm let the symbol x denote a typical time-space point, u the medinm 4- velocity, T the 
absolnte temperature, K = (iL°, K) the gluon 4-momentnm, 16 the spin-colonr degen¬ 
eracy factor, \g the glnon fngacity, / the one-body Bose-Einstein distribntion fnnction, 
and Ug the glnon nnmber density given mathematically by 

/ = = ( 6 ) 

It may be noted that generally < 1 for chemically nnequilibrated plasma and K ■ u 
is tedious in the hreball frame. 


Comoving Frame: If k = {k^, k) is the glnon 4-momentnm in the local rest frame of 
the plasma then by Lorentz transformations 

K ■u = k° ] K^k^ = -f{l + v-k) (7) 


where 7 , v refer to the flnid motion of (1) and /c is a nnit vector along k. Snbstitntion 
into ( 6 ) simplihes the distribntion fnnction as 


/ 


A. 

- Ag 


00 


n=l 


^-nkVT 


( 8 ) 


Also, npon taking the polar axis for integration along v and writing v ■ k = |n| cosOkv 
our number density rednces to 

d^k 


, , f d-k K-^ 

Ualx) = 16 / -Q-rTT/ 

^ J ( 27 rf k^ 


(9) 
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4 

Tl^ JO 


16 oo \n 

= 

'TT-^ • ^ rnO 


TT^ 


n=l 


rr 



d cos 9kv 



d(j)kv 7 (1 + 1^1 cosOkv) E ^9^ 

n=l 


( 10 ) 

( 11 ) 


Remarks : This result is new and shows in a compact manner how the number den¬ 
sity depends upon 7 , T, and \g (although these dependences are generally interwoven 
through the evolution equations). It may be stressed that inspite of the occurrence of 
a V ■ k term in (7) our hnal Ug depends on through 7 . Our (11) generalizes an ear¬ 
lier expression obtained by Xu et al [ 8 ] who ignored the flow velocity v completely and 
worked with a factorized form of the Bose-Einstein distribution containing Xg only in 
the numerator. 


Numerical Estimate at tj . The behaviour of (11) can be easily studied at the instant 
when the thermalized hreball was formed in high energy heavy ion collision. The initial 
conditions predicted by HIJING Monte Carlo simulations are summarized in Table 1. 


Table 1: Initial values for the time, temperature, fugacities etc. at RHIC(l), LHC(l) 
only [15] 



T(GeV) 

U = To (fm) 

A. 

A. 



RHIC(l) 

0.55 

0.70 

0.05 

0.008 

1.78 

4.05 

LHC(l) 

0.82 

0.5 

0.124 

0.02 

14.73 

33.63 


There the gluon densities computed via (11) in the nonrelativistic (n ~ 0.1c) and 
nltrarelativistic (n ~ 0.9c) regions are also listed showing a relative enhancement by the 
factor 


(1-0.01)^/^ 

( 1 - 0 . 81 )^/^ 


2.28 


( 12 ) 


Temporal Evolution; It is very tedious to employ the exact distribution function 
/ of ( 6 ) for determining how the fngacities and temperature evolve with the proper 
time r of the medium. Hence we shall directly borrow the master rate equations from 
existing literatnre based upon approximately factorized form of /. It is plausible to 
assume that such an approximation will not markedly affect onr hnal conclusions since 
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Ag < 1. For longitudinal expansion parametrized by (4) the relevant ordinary differential 
equations [10] are known to be 


+ 3— H— — -Rs (1 — \g) — 2i?2 

Xg IT 


Xq T T 


Art Art 


bl V 



T'^t = const 



(13) 


Here Xq(Xg) is the quark (gluon) fugacity, Nf the number of flavours, and remaining 
symbols are defined by 


R 2 — O.Snglvagg - ,gg), — 0.5ng{V(Jgg - ^ggg) 

ai = 16C(3)/7r^ 02 = SttVIS 

61 = 9C(3)iV//7r^ 62 = 77r^Nf/20 (14) 


Next, for transverse expansion parametrized by (5) the appropriate partial differential 
equations [12] read 


+ r-^dr{rT°^) + r“^(T°° + P) = 0 


(15) 


and 


d^T 


+ r-^dr 


^(yUU ^ p -^^2 p 


(16) 


where 


T°o = (e + P)m°m° - P 


(17) 


Their solutions on the computer yield the functions T{x), Xg{x) and hence ng{x) 
subject to the stated initial conditions. 


IV THERMALLY-AVERAGED RATE 

Preliminaries : Next, we turn to the question of applying statistical mechanics to 
gluonic break-up of the J/ijj moving inside an expanding parton plasma. In the fireball 
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frame consider a ip meson of mass four momentum three velocity 

and dilation factor 7 ^ defined by 

1 /2 

= Py/pJ , 77 = “ "4) 

The invariant quantum mechanical dissociation rate T for g — tp collision may be written 
compactly as 


r 


(19) 


where Urei is the relative flux and a the cross section measured in any chosen frame. Its 
thermal average over gluon momentum in the fireball frame reads 


(r(a;)) 


16 r d^K ^ 
ng{x)J [27rf 


( 20 ) 


Here the tilde implies that the gluon is hard enough to break J/ip and the gluonic 
distribution function / is evaluated at the location of the tp meson viz. at the time-space 
point 


X = 


( 21 ) 


The highly nontrivial integral (20) is best handled in the tp meson rest frame by gener¬ 
ating several useful kinematic informations as follows. 


Kinematics in .J/ip Rest Frame: Let q = q) be the gluon 4 momentum measured 
in Ip meson rest frame. Since the relative flux becomes = c = 1 hence our invariant 
T reduces to the QCD [16] based cross section 


T = = 


= - > 1 ; 


B{Q^ - 

3 V 3 y mde^mcp^^ 


( 22 ) 


here is the J/ip binding energy and rOc the charmed quark mass. This possesses 
a sharp peak at the gluon energy 

,0 _ 1!^ _ 0.92 GeV; Q» = “ (23) 

Also, the energy variable for the massless gluon transforms via 


= 77 ( 9 ° + A/, ■ 9 ) = 779 ° (1 + Wp cos 6 * 5 ^) 


(24) 
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with 9q^ being the angle between q and unit vectors. Furthermore, the fluid 4- velocity 
w = {w^, w) seen in ip rest frame will be explicitly given by the Lorentz transformations 


= 7 ^ M = 7^7 [l-v-vp] 

(25) 

w = 

u-{u-vp)vp + 7 ^ {u ■ vp) - u°\vp\ vp 


= 7 [h'-77W + (7y - l)(h'-w)w] 

(26) 

Finally, the scalar product 

K 

■ u = q ■ w = q^w^ — q^\w\ cos9qyj 

(27) 


where dqw is the angle between q and w. We thus have all the ingredients needed to 
calculate the thermally-averaged rate of (20). 

Evaluation of (r(a;)): Upon taking the polar axis for ^integration along denoting 
the solid angle element by dQq^, and expanding / in a power series we can write 


{T{x)) = -( 
Ur, J 


cAtt 


d^q \n -nK-u/T 

[271) q 

(28) 

OO 

Pcos9qp)aKestT.>^>~"'''"^^ 

(29) 


n=l 


TT'^Ug Je^ Jo 

This integral is performed in the Appendix based on the following nomenclature. 
Symbols/Notations: Using the nomenclature in the fireball frame we define 

{ 9 v, <Pv) = polar angles of v ; { 9 v^,<PvJ = polar angles of v^p 
F = v- vpj = \v\ sin 9^ sin 9^^ cos(0^ - (pv^) + cos 9y cos 9y^ 

y = 7ylwl - (7y - 1 )^ 

= 77y(l - y\vpj\) ; w = 7(F- Yvpj) 


w = w 


1 27/2 


= 7 


\v\" + Y^ 


2YF 


1/2 


9p,w = angle between vp, and w 
1 1 

cos6 '^^ = = t-^(E - Y) 


w 


w 


(30) 


(5° = g°/e^ ; = nep.w'^/T 

Ffl Tl(lqp\w\lT , Pyi DyiCp 
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Also will be needed 



= Cn± Dn = ± |w|)/T 

Io{Pn) = An\l{pn)/Pn 
h{pn) = cosh(p„)/p„ - sinh(p„)/p^ 
nK ■ u/T = nq ■ w/T = CnQ^ — PnCosOqyj (31) 


Result of Integration; From Appendix (A4) we find 


(r(a;)) 


847p 

TT^ng 


OO n QQ 

dgVVRest(g°)e-^"^° 

n=l ^ 


X 


-^o(Pn) A A(Pn)|^p| COS 


(32) 


By decomposing the hyperbolics Jq and Ji into exponentials this result can also be 
expressed as 


(r(a;)) 




TT^ tc n 


3 n=l 


n Jl 




X 


(^1 + Vp- t&(l - 1/pn) ^ e ^"*5° 
t(;(l + l/p„) ^ 


(33) 


Remarks : The expressions (32, 33) are original and demonstrate how the mean dis¬ 
sociation rate (r(a;)) depends on the hydrodynamic flow through |t(J| (or w^) as well as 
the angle 9pw in the notation of (30). The structure of exp(±iA„g°) tells that numer¬ 
ical treatment of (32) is convenient if 0 < |t(;| < T/ep while that of (33) is suitable if 
T/ep < |tc| < OO. From the analytical viewpoint it is much more advisable to work with 
the modihed rate 

(f (i)) = n,(i){r(i)) (34) 

This is because (F) is devoid of any Ug factor appearing in the denominator of (32). 
Furthermore, (F) will be seen to enter directly the survival chance e~^ in (51) later. 

Analytical estimate The physical interpretation of our graphs below will be facilitated 
by deriving a rough approximation for (F) as follows. In Eq.(32) we retain only the n = 1 
term of the summation, and approximate the integrand by its peak value at Qj) (cf.23) 
so that the desired estimate of (34) becomes 

{f(x)) « ^ I" dQ«Q«\_, H 

oc Xg'jpH (35) 
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Here the entire dependence on the flow velocity w is contained in the function 

H = loiDiQl) + hiDiQl) \ \ cosO^^ (36) 

with the coefficients Ci, Di and Af read-off from (31). We are now ready to discuss 
some consequences of (35) in three cases viz. static medium in the fireball frame, no flow 
in the J/'ip rest frame, and ultrarelativistic flow in either frame. 

Static Medium in Fireball Frame : Remembering the notations (30, 31) consider 
the hypothetical case 

n = 0; 7=1; F = 0; F = 7^|h'^| 

= 7 ^ ; w = -7^1^ ; Cn = ne^ 7 ^/T, cos 6*^^ = -1 

Dn = aiwl ; = Cn{l ± |hVI) (37) 

This is precisely the case treated in our earlier paper [9] where the g — ip relative flux 
was correctly taken in the ip meson rest frame. Eqs.(32, 33) also improve the work of Xu 
et al [8] who had written the g — ip relative flux in the fireball frame. Figs.l, 2 borrowed 
from Ref. [9] display the temperature and transverse momentum dependence of the usual 
rate (T) based on (32, 33) relevant to the LHC(l) initial fugacity A^j. For the sake of 
comparison the corresponding curves of the modified rate (F) using (34) are drawn in 
Figs.3, 4. Due to the assumed absence of flow there is no inhomogeneity with respect 
to X and the .J/ip formation time is also ignored at this stage. Various features of the 
indicated diagrams are interpreted in the next paragraph. 

Interpretation i) Following the arguments given by Xu et al [8] it is known that the 
peak of at gives a rich structure to the usual rate (F) in Figs 1, 2. ii) However, 
the modified rate (F) in Figs.3, 4 becomes structureless , i.e., monotonic. This is because 
the gluon number density Ug in (11) contains a crucial factor which changes from a 
low value viz. 0.008 at T = 0.2 GeV to a high value viz. 1.0 at T = 1 GeV. Such a 
coefficient matters a lot in the conversion of (F) to (F) via (34). iii) At fixed pt the 
steady increase of (F) with T in Fig. 3 is caused by the growing exp {—AfQ^) factors 
of the estimate (35, 36). iv) At fixed T the monotonic decrease of (F) with px in Fig. 
4 has a very interesting explanation. For the situation (37) under study w = —is 
antiparallel to so that cosO^w = — 1- Hence partial wave terms Jq and A of H (36) 
interfere destructively in Fig. 4 making (F) small as | A/, | grows. 

No Flow in J/ip Rest Frame: Next, consider a configuration in the fireball frame 
such that the 3 velocities of the plasma and ip meson coincide at some x. This is 
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LHC(l) 



Figure 1; The thermal-averaged gluon--^ dissociation cross section (F) = as a 

function of temperatnre at different transverse momenta px as in [9, Eq.(19)], i.e., in 
the absence of longitudinal flow. The initial glnon fugacity is given in Table 1 at LHC 
energy. 

possible if, for example, both the fluid and ip are moving in the transverse direction. 
Then in (30, 31) we put 

v = v^] F = \v^\ ] = 1 ; uJ = 0 

Cn = ne^/T ; Pn = 0 ; /o(Pn) = 1 ; h{Pn) = 0 (38) 

Inserting this information into the exact formulae (32, 34) and attaching a suffix 0 we 
hnd 

(f (i)>„ = f; a; IP dQO (39) 
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Figure 2: The variation of (F) = with transverse momentum at different temper¬ 

atures [9] as in Fig.l in absence of flow prohle. 

Its variations with T and pt are displayed in Figs.5, 6 relevant to the LHC(l) initial 
fugacity. The crude proportionality (35) now reduces to 

(f(a;))o oc (40) 

which is utilized below to explain some features of the graphs. 

Interpretation i) At hxed pt our (F(a;))Q of Fig.5 rises monotonically with T in analogy 
with the earlier Fig.3. This is caused by the factor present in the estimate (40). 

ii) At hxed T our (F(a;))Q of Fig .6 grows steadily with pt in contrast to the earlier Fig.4. 
Such a behaviour is due to the coefficient 7 ^ occurring in (40). iii) The curves of Figs.5, 6 
are consistently higher than those of Figs.3, 4. The reason is that there is no Ii{pn) term 
present in (39) to interfere destructively with the Io{pn) term in view of the restriction 
(38). 
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LHC(l) 



Figure 3: The variation of the modified rate (f) = ng{r) with temperature at different 
values of transverse momentum corresponding to (F) of Fig.l is depicted here. 

Ultrarelativistic Flow in Either Frame : Finally, suppose at a point x in the fireball 
frame the medium is flowing ultrarelativistically (| F |—;► 1) and the J/fi is moving slowly 
(|F/;| < 1/10, say). In the ip meson rest frame the Lorentz transformation (26) shows 
that w almost equals yu so that cosO^w -v. [Hence, in the case of pure transverse 
expansion of the plasma cos can even become +1 implying a constructive interference 
between the Iq and terms of (36); this possibility will, however, be discussed in a future 
communication. At present, we illustrate the case of both the J/ip and plasma moving 
ultrarelativistically (in the transverse and longitudinal directions, respectively) subject 
to the following kinematic conditions: 

v = 0.9 e^; 7 = 2.3; F = 0; Y = 

\ 

= 77p >> 1; |F1| ~ 77p - 7 :;- 

2777 
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LHC(l) 



Figure 4: The variation of the modified rate (f) = ng{r) with transverse momentum at 
different values of temperatures corresponding to (F) of Fig.2 is shown here in absence 
of flow. 


cos6^^, =^Di = e^|w|/T ^ >> 1 


w 


7 


T 






2 ^ 77 ^ 

Then /o(pi) —^ exp(pi)/2pi so that our rough estimate (35, 36) become 


(41) 


IV( W ^9^ 

(f (x)) oc —— exp 

7 


2 ^ 77 ^/ 


7 


(42) 


This information will be utilized below for explaining the main features of the graphs. 


Interpretation i) For hxed pti v the exponential in (42) tend to 0 as T —0 and 
tends to 1 as T —> cx). Therefore, the growing trend of (F(a;)) with T in Fig. 7 is 
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Figure 5: The variation of the modified rate (f (a;))Q from Eq.(39) with temperature for 
different pfis, when the 3 fluid velocity w in the ip meson rest frame is zero. 

understable, ii) At hxed T, v the rich behaviour of (r(a;)) with pt in Fig. 8 arises from 
a sensitive competition between the bracketed factors of (42). iii) More precisely, at lower 
temperatures T, /eO.4 GeV the exponential factor in (42) increases dominantly with pt 
causing (F(a;)) to grow, iv) But at higher temperatures T > 0.8 geV the third bracket 
in (42) decreases prominantly with px causing (F) to drop, v) Comparison of Figs.(3,4) 
with Figs. (7,8) tells that the pt- dependence of (F) (unlike its T dependence) is quite 
sensitive to the nonrelativistic/ultrarelativistic nature of flow in the hreball frame. 
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p_T (GeV) 


Figure 6: The variation of the modified rate (f (a;))Q from Eq. (39) with pfis for different 
values of T’s when the 3 fluid velocity w in the ip meson rest frame is zero. 

V J/p SURVIVAL PROBABILITY 


Next, we ask the important question as to how the explicit collective flow prohle affects 
the birth and death scenarios of tp mesons within a quark gluon plasma? Below we 
consider only longitudinal expansion; the case of transverse flow will be dealt with in a 
future communication. 

Pure Longitudinal Expansion; Suppose at general instant t in the hreball frame 
the plasma is contained inside a cylinder of radius R, cross section A = nR^, length L, 
and volume V = AL. Keeping the origin at its centre and assuming constant speed of 
longitudinal expansion we have 

R = Ri, A = Ai, L = Lit/U, V = ipt/ti (43) 
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Figure 7: The variation of the modified rate (f (x)) using Eqs.(32, 34) as a function of 
temperature at different transverse momenta for the ultrarelativistic longitudinal flow 
velocity v = 0.9 c. 


where the suffix i lebels initial values. The ansatz (4) of the velocity prohle reads 


V = zezjt ; — T/2 < 2 : < +L/2 


(44) 


subject to the restriction that the speed at the edges ±L/2 must be less than c, i.e., 


2t 2ti 


1 


(45) 


Production configuration of J/fi: Employing cylindrical coordinates let a typical 
meson be created at the instant tj, having position and with transverse velocity 
(in the mid-rapidity region) such that 

ti = ti + ; Tj, = (r^, 4) = (4’ 4) 
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p_T (GeV) 

Figure 8: The variation of the modified rate (r(a;)) using Eqs.(32, 34) as a function 
of pt at different values of temperatures as in Fig.7 corresponding to ultrarelativistic 
longitudinal flow velocity v = 0.9 c. 

vlp = v^ = 0 ) = ( ut , 0 , 0 ) 

F = v- v^ = {zejt) ■ = 0 (46) 

where Tp ~ 0.89 fm/c is the formation time of the bound state in the cc barycentric 
frame. Of course, the concept of formation time was not utilized while drawing the 
graphs 1-8 in Sec.IV. 

Kinematics of J/f) Trajectory; The position 3 vector at general instant t becomes 

= [f^, z^] = [f^ + {t- ti)v^, 4] (47) 

The transverse trajectory will hit the cylinder of radius R = Rj after a time interval tpi 
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by covering a distance dm such that 

= Ri ] tRi = dm/vT 

dm = -r^ cos + \jR]-rlp^ sin^ (48) 

In the hreball frame the full temporal range of our interest is obviously 

ti <t<tii ] tji = mm{ti + tm, bife) (49) 

This prescription was also utilized in [9] and it improves the work of [ 8 ] by incorporating 
the formation time. 

Formulation of S{pt)'- Let us return to the modihed dissociation rate derived in 
Sec.IV. The value F = 0 (cf.(46)) greatly simplihes the kinematics of (30). Furthermore, 
since the time-space point x in (21) is to be chosen on the trajectory itself hence 
the notation (r(a;)) of (34) is equivalent to 

{fl«l> = {f(«,p„4)) (50) 

This depends parametrically on in view of the longitudinal flow prohle (44), and the 
upper limit tu {cf. (49)) of the time variable depends on the production conhguration 
rip, (plp. Then by the law of radioactive decay without recombination the effective survival 
chance of a chosen ip meson will be given by the exponential e~'^ with 

W = T'dt ng[t] {f[t]) (51) 

J tj 

Upon averaging over the production conhguration of the ip'’s we arrive at the dual 
expression for the net survival probability 

S{pt) = [ d%{Rj - rlp‘')e-^/ [ d^4(i??-rf) 

JVj JVj 

d^xlp = drip rip d(plp dzlp (52) 

Remarks : This is a new result showing that due to the inhomogeneity introduced by 
the ensuing longitudinal expansion the spatial integration in (52) must extend over the 
volume Vi = Viti/ti available at the instant of creation. Hence the entry concerning the 
initial length Lj and radius Ri of the plasma becomes essential in Table 2. However, in 
the absence of collective how the integral needs to run only over the cross sectional area 
H/ of the hreball as was done in Refs [9, 8 ]. 
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Table 2: Colliding nuclei, collision energy, initial length of cylindrical QGP, and its 
radius at RHIC(l) and LHC(l). The length is assumed to lie in the range 0.1 < Li <1 
fm since the sea-quarks of the nucleon are spread over a distance of order A 1 fm. 
[We do not use the Li corresponding to directly Lorentz-contracted, disc-shaped nuclei 
to avoid too large or too small values.] 


Nuclei Energy (yi) Li Ri 

(GeV/nucleon) (fm) (fm) 


RHIG(l) 


200 

0.1 <Li<l 

6.98 

LHG(l) 

p^208 

5000 

0.1 <Li<l 

7.01 


Numerical Illustration : As pointed out in Ref. [9] the precise value of the ip meson 
formation time is ambiguous; however we adopt ~ 0.89 fm/c as a fair repre¬ 
sentative. For chosen creation configuration of the tp meson the function W was first 
computed from (51) and then the survival probability was numerically evaluated us¬ 
ing (52). Figs 9 and 10 show the dependence of S{pt) on the transverse momentum 
corresponding to the LHC(l) and RHIC(l) initial conditions, respectively (cf. Tables 
1, 2). The dotted (Lj = 0.1 fm) and dashed (L* = 1 fm) curves are computed in the 
presence of longitudinal flow while the solid curve is borrowed from [9] in the absenee 
of hydrodynamic flow prohle. The main features of these graphs may be explained as 
follows. 

Interpretation i) In the absence of flow, i.e., for n = 0, the S{pt) solid lines in Figs.9, 
10 rise steadily with p^- This is because (T) of Fig.4, and hence the W integral of 
(51), diminish with p^- h) In the presence of flow with small initial length Li ~ 0.1 
fm the S{pt) dotted curves in Figs.9, 10 are slightly above the solid line. To explain 
this we note that the flow prohle v = z/t remains nonrelativistic over short lengths and 
the destructive interference between Jq and R terms in (36) becomes slightly stronger 
(compared to the n = 0 case). The consequent drop of (T) or W pushes S{pt) upwards, 
hi) In the presence of how with large initial length Lj = 1 fm the S{pt) dashed curve 
lies below the solid line in the LHG case (cf. Fig.9) but lies above the solid line in RHIG 
case (cf. Fig. 10). iv) To explain this peculiar contrast between Figs. 9 and 10 we 
note that as the expanding cylinder becomes much bigger than 1 fm the how prohle 
V = z/t may become relativistic over its substantial length so that the graphs of Fig .8 
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can be employed. At higher temperatures (0.8 < T < 0.6, say) at LHC the (L) curves 
are known to drop with increasing But at lower temperatures (0.6 < T < .2), say) 
at RHIC the (L) curves rise with px in Fig.8 whereby a contrast occurs, v) Of course, 
the relative shift in S{pt) due to longitudinal flow is only about 10% to 15% even at 
Pt = 10 GeV. 


0.8 

0.7 

^ 0.6 

CL. 

^ 0.5 
0.4 

0.3 

0 2.5 5 7.5 10 

p_T (GeV) 

Figure 9: The survival probability of J/ijj in an equilibrating parton plasma at LHC 
energy with initial conditions given in Table 1. The solid curve is the result of Ref. [9], 
i.e., in the absence of flow while the dotted and dashed curves represent the survival of 
J/ijj when the plasma is undergoing longitudinal expansion with the initial values of the 
length of the cylinder Lj = 0.1 fm and 1 fm, respectively. 


LHC(l) 
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RHIC(l) 



Figure 10: Same as Fig.9 but computed for RHIC initial conditions. 

VI Conclusions 

(a) In this paper we have extended our earlier work [9] by explicitly including the effect 
of hydrodynamic expansion prohle on the gluonic dissociation of J/ip^s created in an 
equilibrating QGP. The treatment of Secs.Ill and IV is very general in the sense that 
both the fluid velocity v and ip momentum are arbitrary. Only in Sec.V we specialize 
to longitndinal flnid expansion and transverse motion of ip. 

(b) Our theoretical formulae for the gluon nnmber density ng(cf.ll), modihed g — 
Ip breakup rate (F) (cf.34), and ip survival probability S (cf.52) are new. These are 
derived by making careful Lorentz transformations between the rest frames of the hreball, 
plasma, and ip meson. 


22 




(c) At specified fugacity \g the effect of the flow velocity v is to increase the number 

density of hard gluons (which are primarily responsible for breaking the as shown 

in the numerical estimate (12). 

(d) Our expression (32, 34) of (r(a;)) contains partial wave contributions called Iq 
and Ji whose mutual interference is controlled by the anisotropic cos factor. This 
signihcantly affects the variation of (T) with T, px and v as depicted in Figs. 3-8. 

(e) The presence of longitudinal expansion pushes up the graph of S{pt) compared 
to the n = 0 case for nonrelativistic flow appropriate to small initial length Lj = O.lfm. 

(f) However, in the case of relativistic longitudinal flow appropriate to larger initial 
length Li = Ifm the downward shift of S{pt) graph at LHC is in sharp contrast to 
the upward shift at RHIC. Such a contrast between the behaviours at LHC and RHIC 
is caused by the different initial temperatures generated therein. The relative shift is, 
however, is only 10% - 15% even at pr = 10 GeV. 

(g) In a future communication we plan to study the detailed effect of transverse 
expansion of the medium on S{pt)- It is expected that there will be a more rapid 
cooling with time and also a possible constructive interference in (35). 
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Appendix 


EVALUATION OF THE INTEGRAL (28) 


Our derivation proceeds through the following steps : 


First Step We recall the symbols (30, 31), work in the ip meson rest frame, and rewrite 
(29) as 


2 fOO ^ r^TT 

{T{x)) = —— / eldQ^Q^ / 7^, 1 + |h)/,| cos 

TT'^Tig Jl ^ Jo L 

CXD 


(A 1) 


n=l 


Second step Next, a partial wave expansion is done of the exponential 

OO 

gP„cos 0 ,„ ^ + l)Pi{cos9g^) 

1=0 

OO I 

= E Mpn)4n E yr’ (a 2 ) 

/=0 m=—l 

Here pn = as before, ji denotes spherical Bessel functions, Ii{pn) = the 

addition theorem has been used for the Legendre polynomial Pj(cos6*g^), and are 
the polar angles of vp with respect to the local flow direction w. 


Third step Next, the relevant angular integral appearing in (Al,2) reads 




dVL, 


qp 


1 + |n^| cosdg^ Yi^*{VLqp) 


^10 + Y -^\vp\6ii 


= 

due to the orthogonality of spherical harmonics. 


-^mO 


(A 3) 


Fourth step Finally, inserting the informations (A2,3) back into the starting expression 
(Al) we obtain 

O p OO OO 

{r(i)> = — / E 

n=l 

OO 

47rE^z(Pn) 

1=0 

which indeed coincides with (32) of the text. 


dio + 






(A 4) 
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